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Abstract 



o 
o 

CN 1 The differential cross section for scattering of a Dirac particle in a black 

(— i , hole background is found. The result is the gravitational analog of the 

Mott formula for scattering in a Coulomb background. The equivalence 
principle is neatly embodied in the cross section, which depends only on 
the incident velocity, and not the particle mass. The low angle limit agrees 
with classical calculations based on the geodesic equation. The calculation 
employs a well-defined iterative scheme which can be extended to higher 
orders. Repeating the calculation in different gauges shows that our result 
^ ' for the cross section is gauge-invariant and highlights the issues involved 

os : 

in setting up a sensible iterative scheme. 

m 
o 

\Q ' PACS numbers: 11.80.-m, 03.80.+r, 04.70.Bw, 04.62. +v 

O 

1 Introduction 

o 

Scattering of a charged fermion in a background Coulomb field is a widely 
^ | studied process which can be analysed perturbatively in quantum field theory. 

To lowest order the scattering process is summarised in the Mott formula for 
the unpolarised differential scattering cross section O, 



da 



Z 2 a 2 



Mott 4p2 w 2 sin 4 ( /2) 



1- v 2 sin 2 ((9/2)). (1) 



Here p is the momentum, v — \p\/E, a is the fine structure constant and 
the source has charge Ze. Curiously, there have only been sporadic attempts to 
repeat the analysis leading to the Mott formula for the case of a black hole [|[ || . 
The problem of scattering by a black hole has certainly been tackled by many 
authors (see for example the books by Futterman, Handler and Matzner j| and 
Chandrasekhar Q, or the recent article by Andersson and Jensen |6|). But 
few authors tackle the problem in perturbation theory and there is a notable 
reluctance to analyse the fermion case, with most work carried out for the case 
of a massless scalar field. 

Unlike scattering in a Coulomb field, black hole scattering is complicated 
by the additional effects of absorption and emission. Absorption is due to the 
singularity in the gravitational field and manifests itself as a lack of Hermiticity 
in the fermion wave equation |9j . Emission is due to the Hawking radiation 
and cannot effectively be treated without some of the apparatus of quantum field 
theory |l(J. Despite these complicating factors, we show here that an iterative 
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scheme can be set up, based on the Dirac equation, which produces a formula 
for the lowest-order scattering cross section with little difficulty. 

In this paper we concentrate on the scattering cross section for a fermion in 
the background field of a spherically-symmetric black hole. We find that the 
unpolarised differential scattering cross section is given by 

i " iSpC + 2 « 2 - 3 » 2 ™ 2 W2) + "* - «* 8i " 2 W2)) - (2 > 

We calculate this result in two quite different gauges, which require going to 
different orders in perturbation theory. The first calculation utilises a series of 
gauge choices for the gravitational fields which convert the Dirac equation into a 
simple, Hamiltonian form |^|, ^. This can be analysed straightforwardly through 
a perturbation expansion. Unlike approaches based on the Schwarzschild coor- 
dinate system, the gauge used here does not contain a singularity at the horizon 
and extends right up to the origin. This means that each step in the perturba- 
tion series can be treated exactly, avoiding the problems encountered by Collins 
et al. H who attempted a Born approximation scheme based on Schwarzschild 
coordinates. The result they obtained was physically unreasonable, though their 
explanation of why this occurred (as the result of wavepacket dispersion) seems 
incorrect in light of the present calculation. 

The choice of gauge employed in our first calculation leads to an unusual 
form of the vertex factor in momentum space, which vanishes when both the 
incoming and outgoing fermions are on-shell. It follows that the first-order con- 
tribution to the scattering cross section in this gauge is identically zero. Since 
the vertex factor turns out to go as the square root of the black hole mass, the 
fact that the process is second-order does make sense. The integrals involved 
in the second-order calculation are all finite and do not require any form of 
regularisation procedure. The result of the calculation (J2J) is the gravitational 
analog of the Mott scattering formula. For low energies the cross section reduces 
to the familiar Rutherford formula. As the energy is increased, relativistic cor- 
rections become more significant. The formula neatly embodies the equivalence 
principle, in that the cross-section depends only on the particle's velocity, and 
not its mass. In the low- angle limit our formula agrees with earlier results for 
the classical cross-section based on the geodesic equation . 

The success of the calculation suggests a more general scheme for tackling 
the scattering problem, and this is explored in the second half of this paper. 
We calculate the cross section in a different gauge, using a first-order scheme, 
and confirm that the same result is obtained. This verifies that the cross-section 
formula is gauge- invariant. The second calculation also highlights the conceptual 
difficultly of carrying through a full programme based on perturbation theory, as 
choice of gauge at the start of the calculation can dictate features as important 
as the order of Feynman diagrams one needs to consider. This appears to be 
the correct explanation for the failure of the Collins et al. Q calculation. Their 
scheme was based on the Schwarzschild gauge, and they only went to first order 
in perturbation theory. But to be consistent we should include terms from the 
Schwarzschild metric to first order in M, which then recovers the correct result. 

Throughout we employ units with c = Ti = 1 , while factors of G are stated ex- 
plicitly. Where appropriate, factors of c and h are also included. The Minkowski 
spacetime metric has signature (1,-1,-1,-1). 
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2 The Dirac Equation 



Our starting point is the Schwarzschild line element in its standard form 

dg2 = / _ ?9K\ d p _ A _ 1 dr 2 -r 2 (d0 2 + sin 2 9 dtf), (3) 

where i is the proper time measured by stationary observers, and r, 9 and 
(j) have their usual meaning. We first transform to a new time coordinate t, 
corresponding to the proper time for particles infalling radially from rest at 
infinity. This coordinate has 

(2GMr) 1 / 2 , 

dt = dt+ - - dr. (4 

r — 2GM y ' 

In terms of this new coordinate the line element (||) becomes 

-2GM\V2 , \2 



ds 2 = dt 2 - 



(dr + (^-) dt) - r 2 (d9 2 + sin 2 9dtf ). (5) 



The coordinate system is valid for < r < oo. We will not consider the effects of 
analytic continuation through the introduction of Kruskal-Szekeres coordinates. 
See ^ ^] for a more detailed discussion of the form of the metric employed here. 

In order to write the Dirac equation in its simplest form it is useful to revert 
to Cartesian coordinates by writing 

x = r sin# cos0 
y = r sin# sin</> 

z = r cos9. (6) 

With the coordinates written = (t, x, y, z), ^ = . . . 3, the line element (||) 
can be written as 

, 2 j v 2GM 2 2 (2GM\ 1/2 i , 

ds = ri^dx^dx dt A x l dx l dt (7) 

r r \ r J 

where i = 1 ... 3, rj^ is the Minkowski metric and repeated indices are summed 
over. To write down the Dirac equation in this background we start with the 
standard Dirac 7- matrices and from these we define the sets and {g^{x)} 

by 



o o (2GM\ 1/2 x l o , , 
9 = 7 . 9=7- —7 (8) 



and 



(2GM\ 



1/2 



x 



9o=l + — li, 9i = li (9) 

\ r J r 

It is straightforward to check that these satisfy 

{g li ,9v) = 2g llv I, {g^g v } = 28 v J (10) 
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where g M „ is the metric specified by (0) and / is the identity matrix. In terms 
of the matrices the Dirac equation is 

ig^V^ = mi>, (11) 

where 

v^ = (^ + |rfs^)v, S^ = l[ 7 „,7/J], (12) 

and the components of the spin connection are found in the standard way 
(see |fLl[ , for example). With our choice of matrices we find that the Dirac 
equation takes the simple form 

o /2GM\ 1/2 / d 3 \ , 
W _ i7 of__j f_ + _j^ = m ^. (13) 

Here $ = 7 M 9 M is the familiar Dirac derivative operator in Minkowski spacetime. 
The gravitational effects are contained in a single interaction term in an anal- 
ogous manner to the Coulomb interaction. The main difference is the presence 
of a radial derivative. This form of the equation is clearly ideal for scattering 
calculations as the interaction can be treated perturbativcly. In the asymptotic 
region the metric coordinates all agree with the standard Minkowski interpre- 
tation, so there is no ambiguity in the meaning of any cross sections computed. 



3 Non-Relativistic Approximation 

Before studying the relativistic cross section it is useful to first consider the non- 



relativistic limit of the Dirac equation. The Hamiltonian form of (13) contains 
a single interaction term 

Hlt/j = ih(2GM/r) 1/2 r- 3/i d r {r 3/4 i;) } (14) 

where dimensional constants have been included. This interaction Hamiltonian, 
which incorporates all (general) relativistic effects, is independent of the particle 
mass and so embodies the equivalence principle. The interaction term is also 
independent of the speed of light, so the non-relativistic approximation of the 
Dirac equation can proceed in the standard manner jjj . There are no spin effects 
to consider, so we arrive at the Schrodinger equation 

h 2 

V 2 V + ih(2GM/r) 1 / 2 r' 3 / 4 d r (r 3 / i tP) = Eijj (15) 

2m 

where we have assumed that we have a stationary state of energy E. To simplify 
this equation we introduce the phase-transformed variable 

* = V cxp (-i(r/a G ) 1/2 ) (16) 

where 
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and is (eight times) the gravitational equivalent of the Bohr radius. The new 
variable W satisfies the simple equation 

_^V 2 $-^* = M. (18) 
2m r y J 

This is precisely the equation we would expect if we used the Newtonian grav- 
itational potential, and the solutions for ^ are Coulomb wave functions. The 
standard arguments about the long range logarithmic phase effects in Coulomb 
wave functions apply equally to the r 1 / 2 behaviour, so the cross section can be 
found in the conventional way ||. So, in the non-relativistic limit, the gravi- 
tational differential scattering cross section reduces to the Rutherford formula. 
We therefore expect that the full, relativistic calculation will give a cross section 
which reduces to the Rutherford formula for small velocities. 



4 Scattering Cross Section 

The Dirac equation (|l3| ) is well-suited to an iterative solution in the standard 
manner. We seek a solution of 

[i^ 2 - B{x 2 ) - ■m]S G (x 2 ,xi) = S 4 (x 2 - x\) (19) 

where 

n /2GM\ 1/2 / d 3 \ 

The iterative solution to this equation is 

S G (xf,Xi) = S F (xf,Xi) + / d A xi S F (xf,x 1 )B(x 1 )S F (xi,Xi) 



+ JJ d x\ d x 2 S F (xf, xi)B(xi)SF{xi 1 X2)B(x 2 )SF{x2 1 Xi) + ■ ■ ■ (21) 

where Sf{x2,X\) is the (position space) Feynman propagator. The interaction 
term B{x) is independent of time so energy is conserved throughout the inter- 
action. Converting to momentum space we define the amplitude 

M =u s ( Pf )Vu r ( Pl ), (22) 

where 

V - B<p flPi ) + J ^B( Pf , k) p*+? +ie B(k,p t ) + (23) 

Here B(pf,p i ) is the spatial Fourier transform of the interaction term, bold 
symbols refer to spatial components only, and for the spinor terms we follow 
the conventions of Mandl and Shaw 13 . In terms of M. the differential cross 
section is given by 

da / m \ 2 . , ,.n . 

3n = yi M l- (24) 
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The Fourier transform of the interaction term is 

B(p 2 , Pl ) = (2GM)^i 1 °J Ae-^ ^ (JL + 1) e*V* (25) 

where bold symbols refer to spatial components only. To evaluate this we first 
write 



B(p 2 ,p 1 ) = (2GM)^(^f(p 1 -p 2 ) ■ 9/(APl P ^ 



A=l 



(26) 



where 



ipx /27r\ 3/2 



IPI 



We therefore find that the momentum space vertex factor is 

flCPa.Pi) = 3^i(GM)^ ^^ /, (28) 

IP2 — Pi I ' 

This vertex factor has the unusual feature of vanishing if the ingoing and out- 
going particles are on-shell, since energy is conserved. It follows that the lowest 
order contribution to the scattering cross section vanishes. This is reassuring, 
as the vertex factor goes as y/M, and we expect the amplitude to go as M to 
recover the Rutherford formula in the low velocity limit. 

Working to the lowest non-zero order in M the transition amplitude becomes 



where 



Now 



M = -97r d GMu fl (p / )7 u /i 7 V(Pi) (29) 



= , d 3 k Pf 2 -k< }£ + m k 2 -p 2 

1 (2tt) 3 \p f -k\V 2 k 2 -m 2 + ie\k- Pl \V 2 ' [ ' 



m 2 = E 2 -k 2 - m 2 = p 2 - fc 2 , (31) 



where E is the particle energy and p 2 = p 2 = p 2 . The pole in the propaga- 
tor is therefore cancelled by the vertex factors, so there is no need for the ie 
prescription. The integral we need to evaluate is therefore 

f d 3 k k 2 p 2 

h = J Wf\p^k\v 2 \k- Pl \v 2 ^ + m ^ 



which is evaluated in Appendix [X]. The result is 

1 

9ir 2 q 

where q — p t — p { . It follows that 



/ i = 7^rr( 2m + 3 ^/+^)-4i?7 ) (33) 



M = -^-^u s { Pf )(2m- 3{f) f + fr) + 8E 1 )u r (p i ) 

= -^^«.(P/)(2£?7 - m)u r ( Pl ), (34) 
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and the differential cross section is given by 
da {2GMm) 2 



u s (p f )(2E 7 -m)u r ( Pi )\ 2 . (35) 



So, despite the complexity of going to second order in the iterative solution, 
the result is quite straightforward. Performing the usual spin sums gives an 
unpolarised cross section of 

do- _ (GM) Tl ^ f + TO )( 2£ ; 7 o _ m )(^. + m )(2£; 7 ° - m)} 



dfl 2q 4 

2 -^Ml{ m 2 {E 2 ^p r p l ) + (2E 2 - m 2 ) 2 + AE 2 Pf - Pl ). W) 



t 

If we let v = \p\/E denote the particle velocity, and 9 the scattering angle, we 
arrive at the simple expression 

m = J§W) {1 + 2v2 - 3 "' e " Am + v ' - v, ' m2m,) - (37) 

The formula has the satisfying property of being independent of the particle 
mass, as one would expect from the equivalence principle. We delay a fuller 
discussion of the properties of this result until after we have established that it 
is gauge- invariant. 

5 The Kerr-Schild Gauge 

The iterative scheme employed here suggests a generalisation to alternative field 
configurations. In effect, what we have done is taken the covariant Dirac equa- 
tion and re-writtcn it in the form 

(ift - m)if> = 0- g^^ip (38) 

and we have interpreted the right-hand side as an interaction term B{x). This 
method will clearly provide a sensible iterative scheme if the right-hand side 
contains a single factor of some power of M . If this is not the case, there can 
be no simple correspondence between the order of the iterative solution, and 
the order of M in the amplitude. The gauge we have exploited to date has 
the feature that B{x) goes as M 1//2 . The obvious question now is whether we 
can do better and find a gauge where the interaction goes as M . This should 
then avoid having to integrate over intermediate momenta in a second order 
diagram. Such a gauge is provided by introducing the Eddington-Finkelstein 
advanced time coordinate, which can be employed to convert the Schwarzschild 
metric to Kerr-Schild form fl4|| , 

ds 2 = r ]llv dx fl dx 1 ' IJvdx^dx", (39) 

r 

where 

l fl = (l,x/r,y/r,z/y). (40) 
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For our g^(x) matrices we choose 



.9° = 7° + — (7°-7r), 9' = f-— -(7°-7r), (41) 



where 



lr = -7i- (42) 
r 

The reciprocal matrices are found to be 

.90=7° 7°-7r, g i =J l 7°-7r. 43 

r r r 

In this gauge the Dirac equation becomes 

._ . iGM „ (d d 1\ 

^ + __ (7 o_ 7r) ^_____j V) = mV)! (44 ) 

which achieves our goal of constructing an interaction term of order M. 
The momentum space representation of the interaction term is now 

B {P2 , Pl ) = GM J d* X e-Vf*l ( 7 » - (-E + ^ + 1) e*V* 

(45) 

The result of this integral is 

2-kGM , n , s AttGM . , 9w , 

S(p 2lPl ) = - -^p-(4E 7 °- jrfi- jfc) - -7^p(P2 -Pi 2 )0fe- A) 

+ ^^(V - Pi 2 )7° - Mfcfe- A)) • (46) 

This form of interaction is certainly not as elegant as our earlier gauge choice, 
but has the advantage that the first-order term in the iterative solution gives 
the O(M) contribution to the amplitude. Since the final two terms in B(p 2 ,p 1 ) 
vanish on-shell, only the first term contributes to the amplitude, and we find 

M = — u s (p f )(2E 7 ° - m)u r ( Pi ), 47 

precisely as obtained earlier. 

This calculation confirms that the differential cross section formula, to order 
M 2 , is gauge invariant. Since a series for the amplitude in this gauge can only 
go in orders of M, this suggests that there can be no M 3 / 2 contribution in the 
earlier gauge setup. That is, the third order term in the expansion of (|2l]) must 
vanish. This is confirmed by (somewhat tedious) calculations. 

6 Discussion 

Our calculations in two distinct gauges have confirmed that, to lowest order in 
the black-hole mass, the unpolarised differential scattering cross section for a 
spin-1/2 particle is given by 

% = + 2 " 2 " 3 " 2em2 ' m + ~ " 4si ° 2( " /2,) - (48) 



As already commented, this formula is independent of the particle mass and 
depends only on the incident velocity. This confirms that the equivalence prin- 
ciple is directly encoded in the Dirac equation, though it remains to be shown 
whether this holds to all orders. The formula also makes it clear that the low 
velocity limit recovers the Rutherford formula. The higher order relativistic cor- 
rections are not obvious, but do agree with the small angle formulae obtained 
by Collins et al. @, who found that as 9 i— > the classical cross section is given 

by 

da AlGM) 2 „ 2 4 . 

Unlike the Collins et al. formula (their equation 12), however, there is no cubic 
term in 9~ 1 in the quantum result, which is an even function of 9. 

The massless limit m h—>- is also well-defined and leads to the simple formula 

da _ (GAf) 2 cos 2 (fl/2) 

dtl sin 4 (0/2) ■ [ ' 

Again, the low-angle limit recovers the classical formula for the bending of light. 
This result also predicts zero amplitude in the backward direction, 9 = tt. Null 
geodesies produce a significant flux in the backward direction, and the fact that 
zero is predicted here is a diffraction effect for neutrinos which goes beyond 
the predictions of geometric optics. A similar prediction of zero back-scattering 
for neutrinos was made in 0]. A more detailed analysis of the cross section 
in the backward direction also reveals a large 'glory' scattering 0, 0]. In the 
geometric optics limit this is attributable to multiple orbits, and in the quantum 
description the glory scattering is described by higher-order terms in GM . To 
describe these effects in the present scheme requires extending to higher order 
in perturbation theory. This is currently under investigation. 

Extending to higher orders also raises the question of the convergence of the 
iterative scheme proposed here. This is not a straightforward issue to address 
as there is no dimensionless coupling constant in the problem. Also, it is not 
clear whether higher-order quantum terms should still be expected to obey the 
equivalence principle. One can easily formulate desirable criteria for conver- 
gence, such as GME < 1 or GMEv < 1, but these are too restrictive, given 
that the low angle formula we arrive at is expected to be valid for all masses 
and velocities. It would appear that the only way to investigate convergence is 
to compute the next order terms in the perturbation series directly. 

This work should also have clarified the importance of working consistently 
to the correct order in M. This is particularly clear in the Schwarzschild gauge, 
where the interaction term contains factors of 1 — (1 — 2GM /r) 1 / 2 . An iterative 
scheme based on this gauge choice should expand out B{x) as a series in M, and 
then keep all of the terms up to the desired order. Such a scheme is workable, 
but has the disadvantage of introducing new vertex terms at each order in the 
series solution. It is straightforward to confirm that such a scheme reproduces 
our result for the fermion cross section, to lowest order. 

We can now explain the failure of the Born approximation discussed by 
Collins et al. 0] for the scalar case. These authors used a similar technique 
of viewing the difference between the true and and flat space metrics as an 
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interaction term, and constructed the amplitude 

T(p 2 , Pl ) = 2GM J <l\r< » I' " (-|^ + £L£*L^ . (51) 

Applying this in the Born approximation as it stands produces an unphysical 
answer, which can be traced to the pole at the horizon. The resolution is that 
the correct scheme involves a series expansion of the amplitude as well so that, 
working to first order in M, we should compute 



T(p 2 , Pl ) = 2GM J d 3 xe~ l( P*-Pi> x + P2 '^f 1 * 



(52) 



Precisely this integral is obtained in both of the the other gauges discussed in 
this paper, up to terms which vanish on shell. Working to first order in M gives 
rise to the differential cross section 

d ° _ (GM) 2 

dn- 4 V *sin\e/2) [1 + v } ' (53) 

appropriate for a scalar field. Again, we see that the equivalence principle is 
obeyed, and the various low-angle and low velocity approximations are retained. 

This work suggests a number of generalisations, the most obvious of which 
is to more general black-hole configurations. In this respect a start for the Kerr 
case has already been made in [^| , where the Kerr solution was formulated in a 
gauge with similar properties to that employed in the first half of this paper. In 
addition, both gauges discussed here look well set up to give a proper, quantum 
description of radiation processes as a particle is accelerated in a gravitational 
field. Classical descriptions of such processes are notoriously tricky and am- 
biguous. A further question is whether the formalism developed here is also 
appropriate for the absorption problem j?], |ie)| . This involves modifying the 
vertex factor in such a way as to explicitly expose the non-Hermiticity due to 
the singularity. We expect to tackle these issues in future papers. 



A Evaluation of I\ 

To evaluate the integral l\ of equation (^2|) we first displace the origin in A;-space 
by (p f +Pi)/2 to get 

We now chose coordinate axis with q and Pf+Pi defining the 3 and 1 directions 
respectively. These vectors are orthogonal as the momenta are on-shell. We 
next introduce the spheroidal coordinates 

fci = a sinhu sinv cos</> 

&2 = a sinhu sim; sin^; (55) 
hi = a coshu cosu 
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where < u < oo, 0<w<7r, < < 27r and a = \q\/2. Exploiting the 
symmetry in k the integral reduces to two terms: 

h = 



1 f°° , r J ■ , • f sinh^M - sm 2 v . , , n s 

l = — 7T~ n / du dv sinhu smw < = = — (»/■+ ®j + 2m) 

2^V i Jo I (sinh 2 u + sin 2 t;) 5 /2 ^ ; 

+ ^-2 J B 7 °)l. (56) 



sinh 2 u sin 2 t> 
(sinh 2 u + sin 2 w) 5 / 2 

These integrals are simple to perform and lead to the result of equation ( ^3|) . 
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